install.packages("ggpubr")

u1 <- c(3.5, 6.5, 4.1, 5.3, 4.9, 4.3, 4.4, 3.5, 5.9, 3.4, 4.6, 4.9, 
                  3.8, 4.3, 5.8, 3.9, 4.6, 3.8, 4.8, 6.9, 4, 6.7, 3.3, 3, 4.2, 
                  3.5, 4.7, 4.3, 4.6, 4.7, 3.2, 5.4, 4.1, 5.4, 4.6, 5.1, 6.4, 
                  3.1, 5.1, 4.9, 3, 3.7, 3, 4.9, 4.7, 4.7, 3.3) 
p1 <- c(16.6, 10.9, 14.2, 17.3, 13.1, 12, 8.3, 11.1, 19.6, 12.6, 14.7, 9.3, 
             12.6, 12.3, 12.7, 11, 12.4, 17, 19, 12.9, 7.8, 9.9, 13.5, 9.8, 21.1, 
             13.6, 11.5, 10.3, 8, 8,7, 18.5, 14.2, 14.7, 11.4, 13.3, 17, 13.3, 
             12.1, 11.1, 15.7, 13.6, 16.2, 16.9, 10.6, 10.3, 9.6, 11.8, 17.3, 11, 9.4)
i1 <- c(595, 462, 509, 485, 475, 469, 392, 488, 499, 533, 338, 480, 
                   350, 411, 296, 318, 459, 797, 151, 394, 243, 511, 176, 660, 
                   374, 237, 503, 207, 313, 323, 326, 360, 214, 428, 652, 367, 
                   353, 202, 525, 426, 423, 683, 246, 262, 477, 271, 314, 393, 408)
b1 <- c(969.1, 617.3, 925.3, 1,139.90, 676, 682.1, 419.3, 725.2, 659.5, 
               944.6, 909, 677.5, 513.2, 602.1, 731.3, 604.2, 723.3, 644.8, 1,049.20,
               512.9, 667, 546.5, 753.9, 583.9, 935.9, 310.7, 534.5, 994.6, 331.4,
               452, 1,069.70, 355.1, 1,212.70, 376.3, 909.8, 960.5, 645.2, 463.2, 
               507.2, 989.8, 338.9, 1,040.90, 917.3, 576.5, 528.9, 417.6, 911.6, 
               634.1, 485.8, 450.5)
o1 <- c(0.00232, 0.0018226, 0.0019391, 0.0019241, 0.0023841, 0.0023772, 
              0.0022469, 0.0026328, 0.0072154, 0.0025399, 0.0025582, 0.0022528,
              0.0017708, 0.0029099, 0.0018304, 0.0017018, 0.0026454, 0.0019551,
              0.0038969, 0.0016662, 0.0026807, 0.0026702, 0.0019569, 0.0016883,
              0.0022432, 0.0017793, 0.0019952, 0.0021086, 0.0020876, 0.0034358,
              0.0022686, 0.0033214, 0.0024197, 0.0018611, 0.0020322, 0.0019719,
              0.0015558, 0.0027916, 0.0023988, 0.0024787, 0.0015866, 0.0024832,
              0.0021339, 0.0017632, 0.0031214, 0.0023125, 0.0016067, 0.0018488, 
              0.0023476, 0.0025673)

u2 <- c(4.3, 4.6, 4.6, 4.9, 5.8, 5.4, 6.4, 4.9)
p2 <- c(12, 14.7, 12.3, 12.7, 17, 13.3, 15.7, 16.9)
i2 <- c(469, 533, 350, 411, 459, 428, 525, 683)
b2 <- c(682.1, 909, 602.1, 731.3, 644.8, 909.8, 989.8, 917.3)
o2 <- c(0.0023772, 0.0025582, 0.0029099, 0.0018304, 0.0019551, 0.0020322, 
        0.0024787, 0.0021339)

u3 <- c(5.3, 4.6, 3.8, 4.3, 5.8, 3, 5.4, 5.1)
p3 <- c(17.3, 12.3, 11, 12.4, 17, 11.5, 13.3, 16.2)
i3 <- c(485, 350, 296, 318, 459, 237, 428, 423)
b3 <- c(1139.9, 602.1, 604.2, 723.3, 644.8, 534.5, 909.8, 1040.9)
o3 <- c(0.0019241, 0.0029099, 0.0017018, 0.0026454, 0.0019551, 0.0019952, 
        0.0020322, 0.0024832)

res1 <- t.test(u1, mu = 4.8, alternative = "two.sided")
res2 <- t.test(p1, mu = 13.6, alternative = "two.sided")
res3 <- t.test(i1, mu = 514, alternative = "two.sided")
res4 <- t.test(b1, mu = 764.1, alternative = "two.sided")
res5 <- t.test(o1, mu = 0.002346, alternative = "two.sided")

cov1 <- t.test(u2, mu = 4.8, alternative = "two.sided")
cov2 <- t.test(p2, mu = 13.6, alternative = "two.sided")
cov3 <- t.test(i2, mu = 514, alternative = "two.sided")
cov4 <- t.test(b2, mu = 764.1, alternative = "two.sided")
cov5 <- t.test(o2, mu = 0.002346, alternative = "two.sided")

bor1 <- t.test(u3, mu = 4.8, alternative = "two.sided")
bor2 <- t.test(p3, mu = 13.6, alternative = "two.sided")
bor3 <- t.test(i3, mu = 514, alternative = "two.sided")
bor4 <- t.test(b3, mu = 764.1, alternative = "two.sided")
bor5 <- t.test(o3, mu = 0.002346, alternative = "two.sided")

res1
res2
res3
res4
res5

cov1
cov2
cov3
cov4
cov5

bor1
bor2
bor3
bor4
bor5

# 2007 t tests

u4 <- c(3.8, 4.4, 5, 4.5, 5.5, 5.6, 5.9, 4.3)
p4 <- c(12, 14.3, 11.9, 12.3, 17.3, 13.1, 15, 16.3)
i4 <- c(465, 563, 350, 426, 512, 442, 524, 669)
b4 <- c(591.4, 950.2, 587.6, 739.4, 652.7, 859.1, 1,025.80, 955.1)
o4 <- c(0.002369, 0.002722, 0.002911, 0.001796, 0.001916, 0.002173, 0.002504, 0.002177)
  
u5 <- c(5.4, 5, 3.8, 4.1, 5.5, 3, 5.6, 4.7)
p5 <- c(17.9, 11.9, 11, 11.2, 17.3, 11.2, 13.1, 15.9)
i5 <- c(502, 350, 291, 312, 512, 243, 442, 424)
b5 <- c(1,131.40, 587.6, 567, 729.9, 652.7, 509.8, 859.1, 1,002.40)
o5 <- c(0.001972, 0.002911, 0.001731, 0.002683, 0.001916, 0.00197, 0.002173, 0.002495)

cov1 <- t.test(u4, mu = 5, alternative = "two.sided")
cov2 <- t.test(p4, mu = 13, alternative = "two.sided")
cov3 <- t.test(i4, mu = 506, alternative = "two.sided")
cov4 <- t.test(b4, mu = 739.1, alternative = "two.sided")
cov5 <- t.test(o4, mu = 0.00242, alternative = "two.sided")

bor1 <- t.test(u5, mu = 5, alternative = "two.sided")
bor2 <- t.test(p5, mu = 13, alternative = "two.sided")
bor3 <- t.test(i5, mu = 506, alternative = "two.sided")
bor4 <- t.test(b5, mu = 739.1, alternative = "two.sided")
bor5 <- t.test(o5, mu = 0.00242, alternative = "two.sided")

cov1
cov2
cov3
cov4
cov5

bor1
bor2
bor3
bor4
bor5

# 2008 t tests

u6 <- c(4.9, 6.2, 6.4, 5.8, 6.6, 6.6, 6.9, 4.9)
p6 <- c(12.9, 16.5, 13.3, 14.4, 18.6, 15.2, 17.1, 17.2)
i6 <- c(467, 540, 351, 442, 492, 449, 519, 639)
b6 <- c(572, 1,038.90, 612.1, 762.8, 675.5, 892.8, 1,026.10, 946)
o6 <- c(0.002384, 0.00272, 0.002899, 0.00186, 0.001885, 0.002249, 0.002534, 0.002174)

u7 <- c(5.2, 6.4, 4.4, 4.4, 6.6, 3.3, 6.6, 6.7)
p7 <- c(18.8, 13.3, 11.8, 13.4, 18.6, 12.3, 11.7, 17.1)
i7 <- c(511, 351, 291, 303, 492, 247, 449, 436)
b7 <- c(1,180.00, 612.1, 547.9, 699.9, 675.5, 492, 892.8, 1,046.00)
o7 <- c(0.002127, 0.002899, 0.001742, 0.00267, 0.001885, 0.001972, 0.002249, 0.002501)

cov1 <- t.test(u6, mu = 6.1, alternative = "two.sided")
cov2 <- t.test(p6, mu = 13.4, alternative = "two.sided")
cov3 <- t.test(i6, mu = 509, alternative = "two.sided")
cov4 <- t.test(b6, mu = 774.5, alternative = "two.sided")
cov5 <- t.test(o6, mu = 0.002463, alternative = "two.sided")

bor1 <- t.test(u7, mu = 6.1, alternative = "two.sided")
bor2 <- t.test(p7, mu = 13.4, alternative = "two.sided")
bor3 <- t.test(i7, mu = 509, alternative = "two.sided")
bor4 <- t.test(b7, mu = 774.5, alternative = "two.sided")
bor5 <- t.test(o7, mu = 0.002463, alternative = "two.sided")

cov1
cov2
cov3
cov4
cov5

bor1
bor2
bor3
bor4
bor5

# 2009 t tests

u8 <- c(7.7, 9.6, 10.1, 10.1, 10.5, 10.2, 11.7, 7.6)
p8 <- c(12.9, 16.5, 13.3, 14.4, 18.6, 15.2, 17.1, 17.2)
i8 <- c(450, 526, 349, 447, 478, 446, 512, 648)
b8 <- c(530.4, 1,000.70, 603, 761.5, 688.5, 902.9, 992.8, 969.4)
o8 <- c(0.00234847, 0.002734045, 0.002858751, 0.001873429, 0.001945444, 0.002162736, 0.002529833, 0.002230806)

u9 <- c(7.3, 10.1, 6, 6.7, 10.5, 4.6, 10.2, 10.5)
p9 <- c(18.8, 13.8, 12.6, 13.6, 19, 12.9, 13, 17.7)
i9 <- c(522, 349, 292, 305, 478, 243, 446, 426)
b9 <- c(1,203.10, 603, 539.4, 690.7, 688.5, 484.3, 902.9, 1,010.90)
o9 <- c(0.00214634, 0.002858751, 0.001766858, 0.002706429, 0.001945444, 0.002009964, 0.002162736, 0.00252493)

cov1 <- t.test(u8, mu = 9.3, alternative = "two.sided")
cov2 <- t.test(p8, mu = 14.6, alternative = "two.sided")
cov3 <- t.test(i8, mu = 509, alternative = "two.sided")
cov4 <- t.test(b8, mu = 731.3, alternative = "two.sided")
cov5 <- t.test(o8, mu = 0.002491831, alternative = "two.sided")

bor1 <- t.test(u9, mu = 9.3, alternative = "two.sided")
bor2 <- t.test(p9, mu = 14.6, alternative = "two.sided")
bor3 <- t.test(i9, mu = 509, alternative = "two.sided")
bor4 <- t.test(b9, mu = 731.3, alternative = "two.sided")
bor5 <- t.test(o9, mu = 0.002491831, alternative = "two.sided")

cov1
cov2
cov3
cov4
cov5

bor1
bor2
bor3
bor4
bor5

# 2010 t-tests

u10 <- c(8.9, 10.2, 10.5, 10.1, 10.2, 10, 11.2, 8.2)
p10 <- c(13.4, 17.9, 13.8, 15.3, 19, 15.8, 18.2, 17.9)
i10 <- c(449, 560, 377, 431, 457, 448, 490, 646)
b10 <- c(520, 998.4, 587.6, 726.7, 698.5, 923.3, 997.9, 909.1)
o10 <- c(0.002295, 0.002845, 0.00278, 0.001806, 0.001964, 0.002102, 0.002453, 0.002195)

u11 <- c(7.9, 10.5, 6.3, 7.2, 10.2, 4.7, 10, 9.8)
p11 <- c(18.8, 13.8, 12.6, 13.6, 19, 12.9, 15.8, 17.7)
i11 <- c(551, 377, 307, 316, 457, 245, 448, 430)
b11 <- c(1,114.90, 587.6, 546.8, 680.1, 698.5, 455.9, 923.3, 1,012.20)
o11 <- c(0.002119, 0.00278, 0.001745, 0.002691, 0.001964, 0.001983, 0.002102, 0.002526)

cov1 <- t.test(u10, mu = 9.4, alternative = "two.sided")
cov2 <- t.test(p10, mu = 15.3, alternative = "two.sided")
cov3 <- t.test(i10, mu = 510, alternative = "two.sided")
cov4 <- t.test(b10, mu = 735.4, alternative = "two.sided")
cov5 <- t.test(o10, mu = 0.002482, alternative = "two.sided")

bor1 <- t.test(u11, mu = 9.4, alternative = "two.sided")
bor2 <- t.test(p11, mu = 15.3, alternative = "two.sided")
bor3 <- t.test(i11, mu = 510, alternative = "two.sided")
bor4 <- t.test(b11, mu = 735.4, alternative = "two.sided")
bor5 <- t.test(o11, mu = 0.002482, alternative = "two.sided")

cov1
cov2
cov3
cov4
cov5

bor1
bor2
bor3
bor4
bor5

# 2011 t tests

u12 <- c(8.3, 9.8, 9.8, 9, 9.5, 8.6, 10.3, 7.9)
p12 <- c(13.5, 19.1, 15, 116, 19.1, 16.4, 18.9, 18.5)
i12 <- c(427, 547, 376, 443, 479, 441, 473, 633)
b12 <- c(503.6, 974.6, 604.1, 775.7, 745, 976.3, 1,002.80, 838.3)
o12 <- c(0.002223862, 0.002780584, 0.002780164, 0.001738504, 0.001906334, 0.002108088, 0.002477067, 0.002172479)

u13 <- c(8, 9.8, 5.9, 6.7, 9.5, 4.4, 8.6, 9.2)
p13 <- c(19.5, 15, 12.8, 13.8, 19.1, 13.1, 16.4, 18.3)
i13 <- c(545, 376, 295, 324, 479, 244, 441, 443)
b13 <- c(1,173.30, 604.1, 568.2, 654.4, 745, 472.9, 976.3, 979.9)
o13 <- c(0.001969669, 0.002780164, 0.001697498, 0.002594651, 0.001906334, 0.001967019, 0.002108088, 0.002515148)

cov1 <- t.test(u12, mu = 8.6, alternative = "two.sided")
cov2 <- t.test(p12, mu = 15.8, alternative = "two.sided")
cov3 <- t.test(i12, mu = 512, alternative = "two.sided")
cov4 <- t.test(b12, mu = 745.7, alternative = "two.sided")
cov5 <- t.test(o12, mu = 0.00245053, alternative = "two.sided")

bor1 <- t.test(u13, mu = 8.6, alternative = "two.sided")
bor2 <- t.test(p13, mu = 15.8, alternative = "two.sided")
bor3 <- t.test(i13, mu = 512, alternative = "two.sided")
bor4 <- t.test(b13, mu = 745.7, alternative = "two.sided")
bor5 <- t.test(o13, mu = 0.00245053, alternative = "two.sided")

cov1
cov2
cov3
cov4
cov5

bor1
bor2
bor3
bor4
bor5

# 2012 t tests

u14 <- c(7.8, 9, 8.9, 8.1, 8.3, 7.4, 9, 6.8)
p14 <- c(13.7, 19.2, 14.7, 15.6, 19.4, 16.3, 18.3, 17.9)
i14 <- c(392, 542, 440, 489, 440, 458, 601)
b14 <- c(504.2, 874.9, 552.2, 728.3, 675.4, 895.9, 954.5, 785.9)
o14 <- c(0.002234977, 0.002768123, 0.002756742, 0.001955268, 0.001877756, 0.001988738, 0.002357981, 0.001947394)

u15 <- c(7.5, 8.9, 5.2, 5.8, 8.3, 4, 7.4, 8.2)
p15 <- c(19.8, 14.7, 12.7, 14, 19.4, 13, 16.3, 17.9)
i15 <- c(494, 282, 325, 489, 247, 440, 438)
b15 <- c(1,081.30, 552.2, 556.1, 650.3, 675.4, 470.8, 895.9, 870.2)
o15 <- c(0.001979119, 0.002756742, 0.001718994, 0.002684402, 0.001877756, 0.00195062, 0.001988738, 0.002543353)

cov1 <- t.test(u14, mu = 7, alternative = "two.sided")
cov2 <- t.test(p14, mu = 16.2, alternative = "two.sided")
cov3 <- t.test(i14, mu = 518, alternative = "two.sided")
cov4 <- t.test(b14, mu = 705.2, alternative = "two.sided")
cov5 <- t.test(o14, mu = 0.002423406, alternative = "two.sided")

bor1 <- t.test(u15, mu = 7, alternative = "two.sided")
bor2 <- t.test(p15, mu = 16.2, alternative = "two.sided")
bor3 <- t.test(i15, mu = 518, alternative = "two.sided")
bor4 <- t.test(b15, mu = 705.2, alternative = "two.sided")
bor5 <- t.test(o15, mu = 0.002423406, alternative = "two.sided")

cov1
cov2
cov3
cov4
cov5

bor1
bor2
bor3
bor4
bor5

# 2013 t tests

u16 <- c(6.8, 8.2, 9.2, 7.5, 8.3, 7.4, 7.6, 6.3)
p16 <- c(13, 19, 14.7, 15.9, 18.8, 16, 18.6, 17.5)
b16 <- c(476.1, 823.2, 452.1, 653, 596.4, 790.2, 857.8, 721.8)
o16 <- c(0.002281689, 0.002730097, 0.002229732, 0.001971721, 0.001935317, 0.002447871, 0.002442896, 0.00205479)

u17 <- c(6.1, 7.1, 4.2, 4.6, 6.5, 3.3, 5.8, 6.5)
p17 <- c(19.7, 14.7, 12.7, 14, 18.8, 13.2, 16, 17.8)
b17 <- c(1,030.10, 452.1, 513.5, 600.4, 596.4, 476.3, 790.2, 785.1)
o17 <- c(0.002010298, 0.002229732, 0.001712051, 0.002643929, 0.001935317, 0.001924275, 0.002447871, 0.002557067)

cov1 <- t.test(u16, mu = 6.5, alternative = "two.sided")
cov2 <- t.test(p16, mu = 15.9, alternative = "two.sided")
cov3 <- t.test(b16, mu = 643, alternative = "two.sided")
cov5 <- t.test(o16, mu = 0.002424661, alternative = "two.sided")

bor1 <- t.test(u17, mu = 6.5, alternative = "two.sided")
bor2 <- t.test(p17, mu = 15.9, alternative = "two.sided")
bor3 <- t.test(b17, mu = 643, alternative = "two.sided")
bor5 <- t.test(o17, mu = 0.002424661, alternative = "two.sided")

cov1
cov2
cov3
cov5

bor1
bor2
bor3
bor5

# 2014 t tests

u18 <- c(5, 7.1, 7.1, 5.9, 6.5, 5.8, 6.4, 5.1)
p18 <- c(12, 18.3, 14.4, 15.2, 19.1, 15.8, 18, 17.2)
b18 <- c(438.2, 756.9, 388.2, 559.3, 526.7, 680, 759.9, 627.8)
o18 <- c(0.002142, 0.002514, 0.003279, 0.001882, 0.001927, 0.002171, 0.002189, 0.002018)

u19 <- c(6.1, 7.1, 4.2, 4.6, 6.5, 3.3, 5.8, 6.5)
p19 <- c(18.9, 14.4, 12.2, 13.6, 19.1, 12.4, 15.8, 18.3)
b19 <- c(835.7, 388.2, 464.4, 545, 526.7, 422.5, 680, 712.2)
o19 <- c(0.001993, 0.003279, 0.001698, 0.002602, 0.001927, 0.001948, 0.002171, 0.002563)

cov1 <- t.test(u18, mu = 6.2, alternative = "two.sided")
cov2 <- t.test(p18, mu = 15.5, alternative = "two.sided")
cov3 <- t.test(b18, mu = 581.5, alternative = "two.sided")
cov5 <- t.test(o18, mu = 0.002474, alternative = "two.sided")

bor1 <- t.test(u19, mu = 6.2, alternative = "two.sided")
bor2 <- t.test(p19, mu = 15.5, alternative = "two.sided")
bor3 <- t.test(b19, mu = 581.5, alternative = "two.sided")
bor5 <- t.test(o19, mu = 0.002474, alternative = "two.sided")

cov1
cov2
cov3
cov5

bor1
bor2
bor3
bor5

# 2015 t tests

u20 <- c(3.9, 5.9, 5.9, 4.8, 5.4, 4.9, 6, 4.5)
p20 <- c(11.5, 17, 13.6, 14.5, 18.5, 14.8, 16.6, 15.9)
b20 <- c(429.8, 649.8, 361.1, 519.8, 503, 596.7, 705.7, 557.2)
o20 <- c(0.002174, 0.002486, 0.002518, 0.001969, 0.001641, 0.001989, 0.002407, 0.002092)

u21 <- c(5.2, 5.9, 3.7, 4.2, 5.4, 3, 4.9, 5.8)
p21 <- c(19.1, 13.6, 12.2, 13, 18.5, 12.6, 14.8, 16.7)
b21 <- c(760.2, 361.1, 476.7, 527.6, 503, 338.7, 596.7, 655.2)
o21 <- c(0.001975, 0.002518, 0.001703, 0.002532, 0.001641, 0.001964, 0.001989, 0.002519)

cov1 <- t.test(u20, mu = 5, alternative = "two.sided")
cov2 <- t.test(p20, mu = 14.8, alternative = "two.sided")
cov3 <- t.test(b20, mu = 559, alternative = "two.sided")
cov5 <- t.test(o20, mu = 0.002442, alternative = "two.sided")

bor1 <- t.test(u21, mu = 5, alternative = "two.sided")
bor2 <- t.test(p21, mu = 14.8, alternative = "two.sided")
bor3 <- t.test(b21, mu = 559, alternative = "two.sided")
bor5 <- t.test(o21, mu = 0.002442, alternative = "two.sided")

cov1
cov2
cov3
cov5

bor1
bor2
bor3
bor5

# 2016 t tests

u22 <- c(3.3, 5.4, 5.8, 4.4, 5.1, 5, 5, 4.6)
p22 <- c(8.8, 14.7, 11.3, 12, 15.1, 13.3, 14.8, 13.7)
b22 <- c(431.4, 614.4, 374.9, 514, 469.6, 575.9, 664.7, 533.8)
o22 <- c(0.001779, 0.002468, 0.002811, 0.001955, 0.002045, 0.002048, 0.002561, 0.001802)

u23 <- c(3.9, 5.8, 3.6, 4, 5.1, 3.1, 5, 4.7)
p23 <- c(15.3, 11.3, 9.3, 12.5, 15.1, 10.4, 13.3, 13)
b23 <- c(795.5, 374.9, 479.5, 494.1, 469.6, 337.9, 575.9, 606.1)
o23 <- c(0.002169, 0.002811, 0.001487, 0.001799, 0.002045, 0.001619, 0.002048, 0.002396)

cov1 <- t.test(u22, mu = 4.6, alternative = "two.sided")
cov2 <- t.test(p22, mu = 14, alternative = "two.sided")
cov3 <- t.test(b22, mu = 520.4, alternative = "two.sided")
cov5 <- t.test(o22, mu = 0.002355177, alternative = "two.sided")

bor1 <- t.test(u23, mu = 4.6, alternative = "two.sided")
bor2 <- t.test(p23, mu = 14, alternative = "two.sided")
bor3 <- t.test(b23, mu = 520.4, alternative = "two.sided")
bor5 <- t.test(o23, mu = 0.002355177, alternative = "two.sided")

cov1
cov2
cov3
cov5

bor1
bor2
bor3
bor5

# 2017 t tests

u24 <- c(2.8, 4.7, 5, 3.5, 4.9, 5, 4.3, 4.3)
p24 <- c(8.8, 14.7, 11.3, 12, 15.1, 13.3, 14.8, 13.7)
b24 <- c(406.9, 530.9, 339.5, 452.1, 453.4, 502.4, 623.1, 473.7)
o24 <- c(0.002149861, 0.002639059, 0.003161029, 0.002195053, 0.001878953, 0.002004836, 0.002328291, 0.002064692)

u25 <- c(3.7, 5, 3.1, 3.6, 4.9, 2.9, 5, 3.7)
p25 <- c(15.3, 11.3, 9.3, 12.5, 15.1, 10.4, 13.3, 13)
b25 <- c(727.7, 339.5, 479.3, 478.2, 453.4, 337.1, 502.4, 576.5)
o25 <- c(0.002003945, 0.003161029, 0.001704334, 0.00248544, 0.001878953, 0.001962837, 0.002004836, 0.002527714)

cov1 <- t.test(u24, mu = 3.8, alternative = "two.sided")
cov2 <- t.test(p24, mu = 13.4, alternative = "two.sided")
cov3 <- t.test(b24, mu = 492, alternative = "two.sided")
cov5 <- t.test(o24, mu = 0.002438039, alternative = "two.sided")

bor1 <- t.test(u25, mu = 3.8, alternative = "two.sided")
bor2 <- t.test(p25, mu = 13.4, alternative = "two.sided")
bor3 <- t.test(b25, mu = 492, alternative = "two.sided")
bor5 <- t.test(o25, mu = 0.002438039, alternative = "two.sided")

cov1
cov2
cov3
cov5

bor1
bor2
bor3
bor5

# 2018 t tests

u26 <- c(3, 4, 4.4, 3.4, 4.2, 4.5, 3.4, 3.9)
p26 <- c(9.2, 13.5, 9.8, 10.8, 14.6, 12.1, 13.9, 12.4)
b26 <- c(375.2, 431.3, 306.7, 377.6, 384.7, 412.2, 579.7, 410.8)
o26 <- c(0.002199, 0.002722, 0.003213, 0.00221, 0.002001, 0.002204, 0.002483, 0.002176)

u27 <- c(3.7, 4.4, 2.6, 3.2, 4.2, 2.9, 4.5, 3.5)
p27 <- c(15, 9.8, 9.2, 8.5, 14.6, 9.6, 12.1, 12.5)
b27 <- c(636.8, 306.7, 352.6, 430.6, 384.7, 271.9, 412.2, 489.4)
o27 <- c(0.00224, 0.003213, 0.00174, 0.00255, 0.002001, 0.001985, 0.002204, 0.002566)

cov1 <- t.test(u26, mu = 3.2, alternative = "two.sided")
cov2 <- t.test(p26, mu = 13.2, alternative = "two.sided")
cov3 <- t.test(b26, mu = 444.9, alternative = "two.sided")
cov5 <- t.test(o26, mu = 0.002412, alternative = "two.sided")

bor1 <- t.test(u27, mu = 3.2, alternative = "two.sided")
bor2 <- t.test(p27, mu = 13.2, alternative = "two.sided")
bor3 <- t.test(b27, mu = 444.9, alternative = "two.sided")
bor5 <- t.test(o27, mu = 0.002412, alternative = "two.sided")

cov1
cov2
cov3
cov5

bor1
bor2
bor3
bor5

# 2019 t tests

u28 <- c(2.7, 3.5, 4, 3.2, 4.1, 4.2, 2.8, 3.5)
p28 <- c(9.2, 13.5, 9.8, 10.8, 14.6, 12.1, 13.9, 12.4)
b28 <- c(348.4, 372.1, 271.7, 323.7, 345.7, 375.5, 533.4, 392.8)
o28 <- c(0.0022, 0.002571, 0.003065, 0.002178, 0.001969, 0.002434, 0.002382, 0.002028)

u29 <- c(3.5, 4, 2.8, 3.2, 4.1, 3, 4.2, 3.4)
p29 <- c(15, 9.8, 9.2, 8.5, 14.6, 9.6, 12.1, 12.5)
b29 <- c(599.6, 271.7, 371.1, 342.7, 345.7, 245.3, 375.5, 437.4)
o29 <- c(0.00237, 0.003065, 0.001753, 0.002575, 0.001969, 0.001963, 0.002434, 0.002599)

cov1 <- t.test(u28, mu = 3.3, alternative = "two.sided")
cov2 <- t.test(p28, mu = 12.9, alternative = "two.sided")
cov3 <- t.test(b28, mu = 430.4, alternative = "two.sided")
cov5 <- t.test(o28, mu = 0.002476, alternative = "two.sided")

bor1 <- t.test(u29, mu = 3.3, alternative = "two.sided")
bor2 <- t.test(p29, mu = 12.9, alternative = "two.sided")
bor3 <- t.test(b29, mu = 430.4, alternative = "two.sided")
bor5 <- t.test(o29, mu = 0.002476, alternative = "two.sided")

cov1
cov2
cov3
cov5

bor1
bor2
bor3
bor5
